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Because of a logarithmic enhancement from soft, collinear 
magnetic gluons, in dense quark matter the gap for a color 
superconducting condensate with spin zero depends upon the 
QCD coupling constant g not as exp( — 1/g 2 ), like in BCS 
theory, but as exp(— 1/g). In weak coupling, the ratio of the 
transition temperature to the spin-zero gap at zero tempera- 
ture is the same as in BCS theory. We classify the gaps with 
spin one, and find that they are of the same order in g as the 
spin-zero gap. 



In cold, dense quark matter, the attractive interaction 
between quarks of different colors generates color super- 
conductivity PHnJ. In this Letter we discuss in what 
aspects color superconductivity differs from the classic 
model of Bardeen, Cooper, and Schrieffer (BCS) [p"l|| , 
and in which aspects it resembles it. 

One way in which color superconductivity differs from 
BCS theory is the dependence of the condensate on the 
coupling constant. In theories with short-ranged inter- 
actions, such as BCS theory, the gap depends upon the 
coupling constant g as the exponential of 1/g 2 . We ar- 
gued previously, though, that static magnetic interac- 
tions are not screened to any finite order in g ^]||. The 
scattering of quarks near the Fermi surface is then loga- 
rithmically enhanced by the emission of collinear, nearly 
static magnetic gluons, and this changes the gap from an 
exponential in 1/g 2 to one in 1/g. 

The explicit value of the gap in weak coupling was 
first computed by Son 0]. Using an elegant renormal- 
ization group argument, he found that there is an insta- 
bility at a scale <po ~ 6oM3 -5 exp(— co/g), where /i is 
the quark-chemical potential, cq = 37r 2 /\/2, and bo is a 
pure number. To explicitly compute the magnitude of 
the spin-zero gap at zero temperature, 4>o, it is necessary 
to solve a gap equation. This was initiated by Son Q. 
In this Letter we first extend Son's analysis to estimate 
the constant bo. To the order in g at which we work, all 
of our results are manifestly gauge invariant. Details will 
be presented elsewhere |T^ |. 

Next, we solve the gap equation at non-zero tem- 
perature, T, and show that the critical temperature 
for the onset of color superconductivity, T c , divided by 
4>a is equal to the value in BCS theory 111]], T c /<po ~ 
0.567+ 0(g). 

Finally, we classify the gaps for massless fermions with 



total spin J = 1. There are two types, longitudinal and 
transverse to the direction of momentum of the quarks in 
the condensate, <fy\ and <p^, respectively. In agreement 
with Son [||, we find that all spin-one gaps are of the 
same order as the spin-zero gap: 4>i/4>o is a pure number 
of order one. 

Our results are of practical importance. Bailin and 
Love assumed in their original analysis that static 
magnetic interactions are screened, so that the gaps are 
BCS-like, and thus tiny, (j>o ~ 1CP 3 /i. Our results are 
only valid perturbatively, but if we extrapolate to strong 
coupling, we find that because the constant bo is huge, the 
gaps can become quite large: as seen in fig. 1, <fio peaks 
at 4>q ~ 0.26 fi, with a big T c ~ 0.15 fi. At AGS energies, 
heavy-ion collisions can probe the region of // ~ 600 MeV 
and T ~ 100 MeV. Therefore, by triggering on collisions 
in which cool, dense nuclear matter is formed, it may be 
possible to observe color superconductivity. 

That J = 1 gaps are not exponentially suppressed is 
important for quark stars. At very high densities, the 
chemical potential of up, down, and strange quarks are 
nearly equal, so the J = color-flavor locked condensate 
is surely favored At intermediate densities, however, 
because of the large strange quark mass, and the require- 
ment of charge neutrality, these chemical potentials will 
differ. This suppresses the formation of J — gaps, 
which are predominantly flavor off-diagonal. The J = 1 
gaps, however, can form between quarks of the same fla- 
vor, and will be significant. 

We follow the notations and conventions of our previ- 
ous work [Q,H • For massless quarks there are four types of 
spin-zero condensates [fj]]4|||: right-handed condensates 
<^V±' an d left-handed condensates (jff^- The superscript 
refers to particles or antiparticles, while the subscript de- 
notes helicity. In perturbation theory, QCD is manifestly 
chirally symmetric, so that the gap equations for q^r and 
are identical order by order in g 2 . Although the mag- 
nitude of the gaps for (j) r and 4>p must then be equal, 
because they are complex numbers, they differ by an ar- 
bitrary phase. This phase represents the spontaneous 
breaking of parity by a spin-zero, color superconducting 
gap in an instanton-free regime p|,|l3| . 

Without loss of generality, then, we can consider only 
the right-handed gaps, denoted as and 4>~~ , and take 
them to be real and positive. Suppressing chiral pro- 
jectors, and the color and flavor indices juj, the gap 
function is 
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HQ) = $ + {Q) A + (q) + 0~(Q) A-(q) . 



(1) 



The condensate is formed from a quark, with four- 
momentum Q = (<7°,q), and a charge conjugate anti- 
quark. A ± (q) = (l±7o7-q)/2 are projectors for energy; 

q = ?4 q 2 = 1- 

Including the gap, from (15) of H the quark propaga- 
tor is 



G{Q) = 



A + (q) 



A-(q) 

<7o - H 1 



(7 ■ Q " M7o), (2) 



where e g is the energy of the quark relative to the Fermi 
surface: 



ef = v / (9TM) 2 +0 ± (Q) 2 



(3) 



The poles with =Fe+ represent quasiparticles and their 
holes, those with =Fe~ quasi-antiparticles and their holes 
H . At the Fermi surface, q = /z, it takes very little energy 
to excite a quasiparticle, e+ = — </> + , and a lot to excite 
a quasi-antiparticle, « — 2/x. 

At one-loop order, from (A35) of || the equation for 
the gap function $(K) is 
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Here Gq(Q) — 1/(7 • Q — ^70) is the bare propagator for 
charge-conjugate quarks |Q. To evaluate the Matsubara 
sum over q° we use spectral representations ]l5[ ]. 

In the gap equation, the gluon propagator A MI/ includes 
the effects of "hard dense loops" (HDL) p5| . The basic 
parameter of the HDL Lagrangian is the gluon "mass", 
m g ; for N c colors and Nf flavors of massless quarks, 
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For the time being we take strict Coulomb gauge for the 
HDL propagator. HDL corrections can be neglected for 
the quark propagator and the quark-gluon vertex, as the 
quark lines are hard, q ~ /i. 

We solve the gap equation by including the effects of 
the superconducting state in the simplest possible way 
for the quark, Eq. (g), and not at all for the gluon. This 
is reasonable in weak coupling, because the scale of the 
condensate, <^>q ~ fx exp(— co/g), is much smaller than 
either /i or m g ~ gfi [fl6|. 

As in strong coupling BCS theory 111]], Q(K) has an 
imaginary part, but for small g this can be neglected 
in QCD p7| . Consequently, the only values of the gap 
functions (/^(Q) which enter into the gap equation are 
those on either the quasiparticle mass shell, (/> + (±e+, q), 
or the quasi-antiparticle mass shell, <fi~(±e~ ,q). 



Gap equations for ± (ej: , k) are derived from 0) via 
projection with A ± (k). As is typical in models of su- 
perconductivity (ll), the dominant terms arise from the 
quasiparticle poles. These correspond physically to scat- 
tering of quarks near the Fermi surface. As this involves 
little energy transfer between the quarks, it suffices to 
use the nearly static limit of the gluon propagator. 

With these approximations, denoting^ = e^, the gap 
equation for </>(fc) = (f> + (ek, k) becomes 
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where bt = 4-\/2 /i/(-\/37r m g ), and b\ = 2 /i/(V3m 9 ). The 
logarithm ~ ln[l/(e 2 — e 2 .)) arises from the cut term in 
the spectral density of a nearly static transverse gluon 
{|,[|. I n t ne gap equation, there are also terms ~ ln(l/g) 
which arise from the non-static transverse gluons and 
from static longitudinal gluons; these produce the con- 
stants b t and b\, respectively. In addition, there are terms 
~ 1 in the gap equation which contribute to the constant 
b' ; we did not compute these terms. In deriving (^) we 
assume that e^, e q < ji, so we introduce a cut-off 6 on the 
q-integration; the final result is independent of 8. 
At T = 0, an approximate solution of (|^) is p^ ] 
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<p(k) = 0o sin(gy k ) , 
9 



(8) 
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where 0o denotes the value of the condensate at the Fermi 
surface, k = fi. (This is similar, but not identical to the 
solution of §,|,0.) As = ln(26/i/0 o ), Eq. (§ requires 
<?2/ M = tt/2, i.e., 



2'v <^1> I 



(9) 



Our results for cq and the prefactor 1/g 5 are in agreement 
with Son |(| . The constants bt and 61 are the same found 
in an independent analysis by Schafer and Wilczek [pj; 
see also Hong et al. [ [To|| . 

In BCS-like theories with zero-ranged interactions, 
such as Nambu-Jona-Lasinio (NJL) models M , all parti- 
cle pairs around the Fermi surface contribute with equal 
weight to build up the BCS-logarithm, so that the gap 
function is constant: ~ g 2 J dq/e q ~ g 2 ln(2§ / (f>a) , with 
solution (f>Q r~j 26 exp(— \/g 2 ). In a model where fermions 
interact with scalar bosons of mass M s ~ g[i [^), scat- 
tering of particle pairs through small angles is favored. 
The collinear singularity is cut off by M s ^ 0, so that 
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logarithmic factors of ~ ln(/i/M s ) ~ ln(l/<?) appear in 
the gap equation. 

In QCD, the scattering of quark pairs through small 
angles is again favored. If the exchanged gluon is elec- 
tric, the collinear singularity is cut off by the Debye mass, 
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This produces ln(l/<?) terms which contribute 



to the prefactor 1/g 5 in 6, Eq. (0). If the exchanged 
gluon is magnetic, the collinear singularity is only cut 
off by the difference in energies between the incoming 
and outgoing pairs. In the gap equation (||), this gener- 
ates the logarithmic enhancement factor ~ ln[l/(e^ — e 2 )]. 
The dependence of the gap function on ej, is then not 
negligible. Quasiparticles with momenta exponentially 
close to the Fermi surface, e q ~ bp exp(— c/g), dominate 
the integral, with a contribution which is enhanced by 
ln(6/i/e 9 ) ~ c/g. The gap function <fi(q) is weighted to- 
wards these pairs, as 4>{q)/(f>Q ~ sin(7rc/2co) ~ 1. For 
quasiparticles which are not exponentially close to the 
Fermi surface, e q ~ fx and c ~ g, the gap function is 
down by <t>{q)/<h ~ g @- 

The temperature dependence of the condensate can be 
computed from Eq. (g) as follows. We assume that the 
temperature T is of the order of the gap at zero temper- 
ature, 4>q. Let us introduce a dimensionless parameter 
k» 1, and divide the integration region into e q > k</>o 
and e q < K(f>Q. Away from the Fermi surface, e q 3> <fio, 
the Fermi-Dirac distribution becomes a Boltzmann dis- 
tribution, so tanh(e g /2T) ~ 1, and thermal effects are 
negligible. Near the Fermi surface, the thermal factor 
tanh(e g /2T) cuts off any singularity, even at the critical 
temperature, T c , when <f>(q) — * 0. Then the gap func- 
tion is the same as (^|) for e/. 3> k<^0; and a constant for 
€k <C K(j)o. Matching the two regions at nfio, and then 
sending k — > oo, we derive the condition 
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where e q = y/(q - y) 2 + 2 (T), with <j>(T) the gap 
at the Fermi surface at a temperature T, and e q — 
\/(<7 — y) 2 + 4>o- This is correct to leading order in 
g. Equation (|o|) implicitly determines the function 
4>(T)/(/)q] it is identical to that obtained in BCS theory 
in weak coupling 0,Q. In particular, the ratio of the 
critical temperature to the zero-temperature gap is the 
same as in BCS, T c /<po — (/2 + 0(g), where the constant 
C = 2eV?r ^ 1.13. (7 ~ 0.577 is the Euler-Mascheroni 
constant.) 

Following the classification of , for massless quarks 
a spin-one condensate has the form 

£ (4 e (Q) • q + 4>t(Q) ■ P(q) ■ 7) V h A e (q) , (11) 

h.e 

where the sum runs over chiralities, h = r,i, and ener- 
gies, e = ±. V r j — (1 ±75)72 is the chiral projector, and 



P(q) = 1 qq a projector onto the subspace orthogonal 
to q. Because a spin-one condensate is a three-vector, 
there are 12 types of condensates, four longitudinal, 
(j)\ e h = </>J[ e • q, and eight transverse, 4>\%, = 4>h e 1 P(<j)- 
This classification is equivalent to that of Q . While the 
spin-zero gaps are symmetric in the simultaneous inter- 
change of color and flavor indices HH^H , the longitudinal 
gaps <fi\ e h are antisymmetric. The transverse gaps fulfill 
a more complicated relationship, ((f)^ e ) T = —<j>^. The 
spin-zero gaps and the longitudinal spin-one gaps do not 
mix quarks of different chirality; the transverse spin-one 
gaps do, and thus break chiral symmetry. The gap equa- 
tions can be constructed as in the spin-zero case 12 . We 



find that both the longitudinal as well as the transverse 
gaps fulfill the same gap equation as the spin-zero gaps, 
with identically the same solution as in @, except that 
the constant analogous to b' (which we do not compute) 
may differ. 

In the static limit, gauge dependent terms in the gluon 
propagator A >il '(P) are ~ p^p^/p 2 (TJ]. These terms 
contribute to the gap equation, but neither to Co, the 
power of g in the prefactor, 6 t , nor b\. They do appear 
to contribute to the undetermined constant b' , but we 
suggest that in the end, b' is gauge invariant. There 
are other effects which contribute to b' fL^f . One-loop 
diagrams with a soft, transverse HDL gluon propagator 
renormalize the quark p8| and gluon wave functions, and 
the quark-gluon vertex. Other contributions arise from 
the influence of the condensate on the gluon propagator 
p6| , and the admixture of quasi-antiparticle modes in the 
quasiparticle gap equation. It is important to calculate 
6g, since its numerical value determines exactly which 
patterns of symmetry breaking are favored. 
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FIG. 1. <f>o/y as function of g, for b' — 1. 

While the results which we have derived are rigorously 
valid only in weak coupling, it is interesting to plot cj>o/fi 
as a function of g, fig. 1. We take Nf — 2 (note that 
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bo ~ l/N^ 2 ). From (^), </>o/m is proportional to 
in fig. 1 we set this undetermined constant equal to 1. 
Equation (0) has the form of a semiclassical tunneling 
probability, including a prefactor from five zero modes. 
Because of the "zero modes" , the gap function peaks at 
a value of 4>q//j. ~ 0.26 when g ~ 4.2. 

Extending the picture of Schafer and Wilczek ||, we 
view quark matter as a color superconducting "liquid" , 
and hadronic matter as a color superconducting "vapor" . 
From |Q there is a first-order phase transition between 
these liquid and vapor phases at /i = // c and T = 0. Then 
perhaps at several times nuclear matter density, the liq- 
uid phase occurs at the maximum of (j)o/ g., and the vapor 
phase at larger g 7 providing a qualitative explanation for 
the smallness of the analogous gaps in hadronic matter 

We conclude by using (|4|) to estimate the validity 
of perturbation theory. Perturbative calculations break 
down when m g ~ g or T. For N c — 3 and Nf — 2, at 
T ^ and g = 0, m g = T when the QCD fine structure 
constant is tiny, a s = g 2 /47r ~ 0.18. In contrast, at [i ^ 
and T = 0, m g = g when a s is much larger, a s ~ 2.4. 
This suggests to us that while perturbation theory is not 
a good approximation for hot quark-gluon matter jn| , it 
may well be a reasonable guide to understanding dense 
quark matter, as long as it is cold, T < 0.3/i. 
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Pisarski and D.H. Rischke 
The interaction between two quarks contains two pieces, 
which are symmetric or antisymmetric in the color in- 
dices of the fundamental representation. The antisym- 
metric representation is attractive to lowest order in g; 
for an SU(N C ) gauge theory, the coefficient in (^) is 
g 2 (N c + 1)/{2N C ) |f. When N c = 3, the antisymmetric 
representation is the color 3 representation, the symmet- 
ric the color 6. Fermi statistics for a J = gap imposes 
constraints which require the number of massless flavors, 
Nf > 2 Q. The form of the quark propagator in ([|) is 
only valid when Nf = 2, and the (SU(S) C , SU r ,e( N f)) 
representation is (3, 1) jS|. When Nf = 3, the (3, 3) rep- 
resentation mixes with the (6,6) M|. This mixing only 
affects the gap equation to higher order ~ which is 
negligible in weak coupling. 

M. Le Bellac, Thermal Field Theory (Cambridge, Cam- 
bridge University Press, 1996). 

Due to infrared singular factors, the effective action for 
the condensate is ~ \D^(f)\ 2 j '4>q, so that the (true) gluon 
mass from color superconductivity is not m 3U per ~ g<f>o, 
as one would naively expect, but much larger, m super ~ 
gg,. (We thank T. Schafer for discussions on this point.) 
To the order at which we work, this is irrelevant for the 
gap equation, because the dominant momenta are S> f^o, 
and on that scale, corrections from the condensate are 
small, ~ g 2 cf>o/q at large q S> 4>o- 

From (0) the imaginary part of (f)(k) arises from the cut in 



the logarithm for e q < tk, Im <fi(k) 

2 



g~ \n{tk/(f)o)4>o- Taking tk ~ bg exp(— c/g), momenta ex- 
ponentially close to the Fermi surface occur when c ~ 1. 
In this region, the imaginary part of the gap function, 
Im0(fc) ~ g(co — c)0o, is down by g relative to the real 
part, Re0(fc) ~ sin(7rc/2co)0o- Away from the Fermi sur- 
face, efe ~ g, so c ~ g, and 0(fc) is strongly damped, with 
the real and imaginary parts of comparable magnitude, 
Re0(fc) ~ Im0(fe) ~ g<f>o. 

T. Holstein, R.E. Norton, and P. Pincus, Phys. Rev. B 
6, 2649 (1973). 
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